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ABSTRACT 

In recent years a plethora of future surveys have been suggested to constrain the nature 
of dark energy. In this paper we adapt a binning approach to the equation of state 
factor "w" and discuss how future weak lensing, galaxy cluster counts, Supernovae 
and baryon acoustic oscillation surveys constrain the equation of state at different 
redshifts. We analyse a few representative future surveys, namely DES, PSl, WFMOS, 
PS4, EUCLID, SNAP and SKA, and perform a principal component analysis for 
the "w" bins. We also employ a prior from Planck cosmic microwave background 
measurements on the remaining cosmological parameters. We study at which redshifts 
a particular survey constrains the equation of state best and how many principal 
components are significantly determined. We then point out which surveys would be 
sufficiently complementary. We find that weak lensing surveys, like EUCLID, would 
constrain the equation of state best and would be able to constrain of the order of 
three significant modes. Baryon acoustic oscillation surveys on the other hand provide 
a unique opportunity to probe the equation of state at relatively high redshifts. 

Key words: cosmology:observations - cosmology:theory ~ dark energy 



1 INTRODUCTION 

One of the most challenging puzzles in modern physics 
is the cause for the observed accelerated expansion of 
the Universe. The renewed interest in accelerating cos- 
mologies was born o ut of Type la Superno v ae observa- 
tions a decade ago ( Perlniutter et al.l 1 19971 : iRiess et al.l 



universes with the observed cosmological constant in a 
range of the multi verses of the meta-stable v acua of the 

string theory landscape (jKachru et al.l |2003|). Th is. in 

connection with the anthropic principle 11 Weinberg 1987 : 
Efstathioul Il995l : IBoussoI l2006l : iPeacockl l2007l : ICline et all 



1 19981 : iPerlmutter et al 



199S). In the meantime the cos- 



mic concordance model ( Ostriker fc Steinhardj 1 1993 ) has 
been confirmed by complementary obs ervations of the Cos - 
mic Microwave Background (CMB) llSpergel et al.l l2003h 
the l a rge scale distribution of galaxies ( Efstathiou et al.l 

|2002|: Tegmark e t al.l 2004 ) a nd c lusters of galax;ies 

I Allen et al. 2004: iRapetti et ahllJOO^ ). In addition obser- 
vations of large scale CMB a nisotropics IjScranton et al.l 
120031 : iPadmanabhan et al.ll2005r ) and Baryo n Acoustic Oscil- 
lations (BAO) in the matter distribution (|Cole et al.l 12003 : 
lEisenstein et al.ll2005l ) also indicate the presence of late time 
acceleration in the Universe. 

The simplest way to model an accelerated expansion 
in the Universe is by including a cosmological constant in 
Einstein's equations of gravity. However, the introduction 
of a cosmological constant requires an extreme fine tuning 
of the initial conditions of our Universe to about 120 orders 
of magnitude. This motivates the search for alternative 
explanations for the speeding up of the expansion rates, 
although there have been recent successes in obtaining 
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2007]), could motivate the observed value of the cosmolog- 
ical constant. Nevertheless, it is still valid to seek other 
alternatives to a cosmological constant. One approach 
is to introduce a dynamical scalar field, usually dubbed 
Quintessence, which is either slowly rolling d own a po- 



tenti al or is tr apped in a false vacuum state (WetterichI 
19881: [Peebles fc Ra tral Il988l: iRatra fc Peebles! Il988|-^ 
Frieman et al. 1993: iFerreira fc Joyce! 1 1998! : IZlatev et al.l 
1999! : lAlbrecht fc Skordis! I2OOO!) . These models can lead 
to cosmic acceleration with less fine tuning. Another 
possibility was that the Universe is permeated by a network 
of cosmic defects, which also acts like a source in the 
energy-mo mentum tensor and can lead to accel erated 
expansion (|Battve et al.! Il999l : ISattve fc Mos3 l2005h . The 
dynamical evolution of the Universe in Quintessence or 
defect based models is governed by the equation of state of 
the dominant components. The equation of state is given 
by the ratio of the pressure to the density of these fluids. 
For a cosmological constant this is -1. However dynamical 
dark energy models can deviate from this value and the 
ratio can also vary with time. Typically the ratio is named 
"w = pdc/pde". The precise measurement of the value of 
"«;" is hence one of the foremost tasks for observational 
cosmology. Any significant measurement of a deviation of 
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"w" from -1, would be a major result. Currently conserva- 
tive constraints on "w" are in the regio n oi w = —1 ± 0.15 
l|Tegmark et al.ll2004l : lAstier et al.]l2006h . 

Up to now, it has been difficult to make any connec- 
tion between suggested Quintessence fields and fundamen- 
tal theories. Another possibility to obtain accelerated ex- 
pansion is the modification of Einstein's gravity on large 
scales. There are two a pproaches along th ese lines. One is 
to modify a proposal bv lStarobinskiil (|l980l ). which includes 
higher order curvature terms in the gravitational action to 
obtain accelerated expansion in the early Universe. This has 
been adapted to explain th e observed accelerated expansion 
l|CarroU et al.ll20o4 . |2005^ . Another approach is to obtain 
effective four dimensional Friedman equations, from a five 
or higher dimensional theory, where our standard model is 
confined to a brane, with only gravity allowe d to leak out 
into the extra dimensio n(s) ( Dvali ot al. 2000; Deffavet et al.l 
I2OO2I : iDvaU fc Turneil[2003h . Although thes e models might 
have structural prob lems such as ghosts (|Kovamal l2005l : 
IPe Felice et al.ll2006h it is a worthwhile exercise to explore 
their cosmological phenomenology. In particular they are dif- 
ferent to Qui ntessence models in the wa y structures grow in 
the Universe (|Sonell2005l : iKovamallJOO^ . Finally it is possi- 
ble, at least in theory, that the observed acceleration is an 
effect, which results from an averaging procedure to describe 
the large scale dynamic s of an intrinsically inhomogeneous 
Universe (|Buchert||2005l '). 

All the models described above have in common that 
they can, but don't have to, deviate from the evolution of 
a so called A-cold dark matter (ACDM) Universe consisting 
of a cosmological constant and dark matter. Effectively the 
background evolution for all these models can be mimicked 
by a suitable chosen additional component in the energy- 
momentum tensor. Hence the quest for "w" goes beyond its 
usual application to Quintessence models, although they are 
the only models where there is actually a physical meaning 
associated with "u)". It is therefore an interesting question 
to ask if "w" is different from -1 the value for a cosmological 
constant. The first attempt to address this question is to look 
for different consfani values of "w" . However, this approach 
is very restrictive given that dark energ y models allow for 
evolving equation of state f actors "w" l|Weller fc Albrechtl 
l200ll . l2002l : lMaor et ai]|2002l '). Hence it is important to allow 
for a time varying equation o f state parameter "w" . A simple 
Taylor expansion in redshift (jWeller fc Albrechtll2002r ) is too 
restrictive if we want to include data sets which are sensitive 
to high redshifts, like the CMB. A more successful approach 
is to parameterise a time varying equation of state with a 
Taylor expansion in the scale factor a, like w = wo+Wa(l— a) 
l|ChevaUier fc Polarskill200ll : lLindeill2003l ). This was also the 
parameterisatio n of choice adapted b y the Dark Energy Task 
Force (DETF) IJAlbrecht et aLlbOOeD . The DETF chose this 
parameterisation to be able to compare different proposed 
surveys on an equal footing. In addition this parameterisa- 
tion allows to some extend the distinction of two subclasses 
of Quin tessence mod els, namely freezing and thawing out 
models (|Linderll2006h . 

In this paper we want to pursue a more general ap- 
proach. In order to figure out the redshift sensitivity to the 
dark energy density, or better its ti me derivative, we will 
bin the equation of state in redshi f t (IHuterer fc StarkmanI 
I2OO3I : ICrittenden fc PogosianI l2005l : iHuterer fc Peirid l2007l : 



lAlbrecht fc BernsteinllJOOTl : ISuUivan et al.ll2007l ). Ahhough, 
there are some ambiguities how to exactly perform the bin- 
ning (Ide Putter fc Linderll2007h . the binning approach is in 
general the most model independent approach to fit for 
the background evolution of the Universe. Note that this 
approac h is inherently related to the weighting funct ion 
method IJDeep Saini et al.ll2003l : ISimpson fc Bridlell2006l ) in 
the limit of an infinite number of bins. 

We are now faced with the difficult task of selecting sur- 
veys to pursue our analysis. We emphasize that our selection 
is subjective and our aim is that we have at least one rep- 
resentative survey for each of the probes we are discussing. 
In this introduction we will mention briefiy the surveys and 
what their technical specification is expected to be, while we 
describe how they exploit particular cosmological probes in 
the sections where we discuss constraint from the different 
probes. 

We start with Stag e II surveys accordin g to the clas- 
sification by the DETF (JAlbrecht et al.l 120061 ). Already on- 
line in its simplest configuration is the PANoramic Survey 
Telescope and Rapid Response System (PAN-Starrs|j. PAN- 
Starrs in it's full configuration will consist of four 1.8 metre 
telescopes, equipped with a wide field camera with a field 
of view of 7 arcminutes and a total of 1.4 Gigapixels. Cur- 
rently only one of the four planned telescopes is deployed. 
We refer to this setup as PAN-Starrsl (PSl) in contrast to 
the fuU setup PAN-Starrs4 (PS4). The PAN-Starrs survey 
is planned to take exposures of ~ 60 seconds in each of the 
grizy filters in the PSl configuration and should go much 
deeper in the PS4 configurations with grizy filters. One of 
the surveys with PSl will be a nearly all sky, 20.000 deg^, 
survey up to a median redshift of z ~ 0.5. The nearly all 
sky survey is h oped to be exploited for large scale struc- 
ture and BAOs (iPeebles fc Yu|[l970l : IHu fc EisensteinlllQQSl : 
iBlake fc Glazebrookll2003l ) experiments as well as weak grav- 
itational lensing observations. The medium deep survey will 
be suitable for cosmological constraints with Type la SNe 
out to redshifts of 2: ~ 0.8. PSl in the DETF category is a 
Stage II survey, while PS4 is Stage HI. 

Another survey already deployed, which is relevant to 
dark energy scien ce, is the South Pole Telescope (SPT) 
JRuhlet alll2004l '). The SPT is a 10 metre telescope with 
a bolometer array in its focal plane, designed to de- 
tect thousa nds of clust ers of galaxies via their Sunyaev- 
Zel'dovich (jSunvaev fc Z eldovich 1972) decrement in the 
sub- millimeter domain over 4.000 deg^. Current plans are 
to operate at 150, 219 and 274 GHz. This survey in addition 
with a redshift survey will allow to constrain dark energy 
with the clustering and redshift distribution of galaxy clus- 
ters. 

We describe now near- future surveys, which are 
currently not deployed. The Dark Ene rg y Sur vey (DES) 
l|The Dark Energv Survev Collaboration! 120051 ) is an 
optical-near infrared survey under construction. The op- 
tical part will go down to 24**^ magnitude in the grizy 
bands on the 4 metre Blanco telescope at Cerro Tololo 
Inter- American Observatory (CTIO). The IR counterpart 
will come from part of the Vista Hemisphere Survey (VHS) 
survey performed by ESA's 4m Visible and Infrared Survey 
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Telescope for Astronomy (VISTA) in the JHK bands to 
~20th magnitudes. The DES camera has a 2.2 deg field 
of view, which allows to obtain imaging and photometric 
informations of millions of galaxies. It is planned to observe 
5.000 deg'^, which overlap with the area of the SPT survey. 
This will allow to constrain dark energy via weak lensing 
tomography, galaxy clustering and galaxy cluster redshift 
distribution. In addition there is a wide-area survey where 
10% of the allocated time will be used to follow Type la 
SNe light-curves. 

The Wide Field M uhiObject Spectrograph (WFMOS) 
iGlazebrook et al.ll2005^ will consist of a new wide-field (1.5 
deg) 4000 fibers spectrograph with a passband of 0.39-1.0 
microns which would be deployed in a 8m class telescope. 
Current suggestions are that there will be 10 low dispersion 
spectrographs with R=1800 in the blue and R=3500 in the 
red. This will enable the measurements of BAOs in 0.5 < 
2 < 1.3 and 2.3 < z < 3.3 using the redshifts of millions of 
galaxies over ~ 2,000deg^ at low z and ~ 500 deg^ at high 
z. 

We will now introduce two proposed satellite missions, 
which might come online in the long term. Both NASA and 
ESA in their "Beyond Einstein" and "Cosmic Vision" stud- 
ies have proposed to pursue missions which are able to probe 
dark energj|f|. One of the contenders for NASA is the Su- 
perNovae Acce leration Probe (SNAP) (|SNAP Collaboration! 
I2OO5I : lAlberTet al. 2005'). SNAP is a 2 metre telescope in 
space with a 0.7 square-degree wide field imager and a 
R ~ 100 spectrograph. Both are sensitive in the 0.4-1.7 /im 
wave-band. SNAP is designed to probe dark energy with a 
SuperNovae and weak lensing survey, where the weak lensing 
takes advantage of the 9 filters with a depth of 26.6AB mag- 
nitude. We will discuss the particulars of the two surveys in 
the relevant sections later in the paper. The EUCLID survey 
is a combination of the former SPACE l|Cimatti et al.|[2()()8l ) 
and DUNE proposals. The Dark UNiverse Explorer (DUNE) 
l|Refregier et al.ll2006r ) is a proposed wide field space imager 
on a 1.2 metre telescope with a 1 deg^ visible/near-IR field 
of view. It is designed to measure cosmic lensing over 20.000 
deg^ of the sky and will exploit this as its main cosmologi- 
cal probe. DUNE is designed to use one broad visible band 
(R-l-I-l-Z) for accurate shape measurement for weak lens- 
ing and Y,J,H in the near infrared to complement optical 
photometry which should be available by 2017 for accurate 
photometric redshifts in the range < z < 2. 

The Square Kilometer Array (SKA) is a radio inter- 
ferometer planned to operate in a large range of frequen- 
cies (60MHz - 35GHz) with a sensitivity of 20,000m^/K in 
the range 0. 5-1.4GHz with a fi eld of view of ~ 200deg^ 
at 500MHz l|CariUi fc Rawlind [20o3 ) . This interferometer 
would allow to locate galaxies in the Universe given their 
21cm line emission and would allow us to pe rform large sur- 
veys for galaxy evolution and Cosmologv iJRawlings et al.l 



I2OO4I : iBlake et al.ll2004l : lAbdalla fc Rawlingsll2005l . l2007l . Ab- 
dalla et al. in preparation). The project is planned to be an 
on-going project which should build up from 1% demonstra- 
tors in the following years to a 10 — 15% core in ~ 5 — 10 
year time to its full completion in ~ 2020. 



This concludes our preview of the surveys being dis- 
cussed in this paper. However we will exploit one more sur- 
vey, not for its ability to constrain dark energy, rather for 
its ability to constrain other cosmological parameters to high 
precision. The forthcoming Planck satellite mission will ob- 
serve the sky in 9 radio wave ba nds in order to measur e 
the anisotropics to in the CMB (|ESA-SCI(2005')1| bOOSl V 
Given that the primary anisotropics are mainly a probe of 
the angular diameter distance to the surface of last scat- 
tering, Planck alone is not a strong probe of dark energy. 
However it puts strong constrain ts on other cosmological 
parameters (|ESA-SCIf2005)lll2005D . We wiU hence use fore- 
casts for the Planck surveyor to put prior constraints on the 
remaining cosmological parameters. 

We start the paper by discussing the method of princi- 
pal components analysis (PCA) in the context of the binning 
of the equation of state in Section [2] In Section |3] we discuss 
constraints arising from the Planck surveyor and introduce 
the covariance matrix we will use in the later sections of 
the paper. In Section |4] we study the principal components 
of Type la Supernovae as cosmological probes. In this sec- 
tion we will also analyse the impact of marginalization of 
the remaining cosmological parameters. In Section [5] we will 
study the PCA in the context of cosmic lensing. Section [S] 
will analyse cluster counts as cosmological probes. In the 
Section [7] we will examine the ability of BAOs to constrain 
dark energy. Before concluding in Section [TD] we will dis- 
cuss how to perform a joint principal component analysis 
between complementary surveys. 



2 PRINCIPAL COMPONENTS OF THE 
EQUATION OF STATE 

We will now introduce the method we use to constrain 
the equation of state of dark energy. As outlined in the 
introduction we will pursue a binning approach to the 
equation of state. Binning in thi s context was first in- 
troduced by iHuterer fc StarkmanI (120031). but has since 
been studied by many aut h ors (ICrittenden fc PogosianI 
2OO5I: iHuterer fc Peiri j l2007l: lAlbrech t fc Bcrn steinI l2o'o7l : 
ISuUivan et al.l l2007l : Ide Putter fc Linder, i2007i '). There are 
different possibilities to bin the equation of state, but the 
one we follow here is given by 



w{z) 



Wi, Zi 2^ ^ 2 ^ .Zi H 2^ 



(1) 



at: 



universe . nasa . gov/home . html 



and 



www.esa.int/esaSC/SEMA7J2IU7E_index_0.html 



where Wi is value of the equation of state of dark energy in a 
given redshift bin z\ ± |Azi. Note that beyond a maximum 
redshift z^ax we assume a constant equation of state factor 
Wh- Although the binning of w in redshift leads to a quasi 
model independent fitting procedure, the increased number 
of parameters in general lead to a better fit but with the 
drawback of greatly increased error bars. Typically we will 
choose the redshift bins in the region of Azi = 0.05, hence 
obtaining tens of new parameters for a given survey. Occam's 
razor tells us that this large immber of parameters do not 
lead to significant improvement of the fit in general. However 
just increasing the bin width or cutting of all surveys at a 
given redshift does not do justice to all surveys we are going 
to discuss. This is because in general the error bars between 
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different w bins are highly correlated. What we want is to 
extract information described in a decorrelated way. This 
can be achieved by diagonalising the correlation matrix of 
the w bins and then expressing the fit in terms of the eigen- 
modes. This is essentially a Principal Component Analysis 
(PC A). By ordering the eigenmodes with respect to the size 
of their corresponding eigenvalues, we build a hierarchy from 
the best constrained modes to the least constrained ones. 

For any given experiment we look at the parameter vec- 
tor, which includes the standard cosmological parameters, 
like the matter contents fim, the baryon contents fib, the 
Hubble constant Ho, the spectral index Ua of primordial per- 
turbations , the amplitude of the primordial power spectrum 
o"8. Note that we restrict our analysis to flat cosmologies. In 
addition, as we shall see later, each experiment has also a few 
nuisance parameters. We introduce the parameter vector 

Q = {6i) = {Q.ui,^h,HQ,...,wi,'W2,'W'i,...,WN,Wh} , (2) 

where we assume A'^ bins in redshift to flt w. In order to 
obtain the correlation matrix for a given survey we need to 
know the likelihood of the data vector x given the param- 
eters G: p(x, 0). From this we can estimate the correlation 
matrix with the help of the Fisher information matrix 



The error bars on the a^ coefficients are then given in terms 



F 



(3) 



with C — —2 In p(x, 0). We will show how to calculate the 
Fisher matrix for the particular surveys in the corresponding 
sections. 

In order to study the ability of a given survey to con- 
strain dark energy in different redshift bins we first have to 
marginalize over the other cosmological parameters and the 
nuisance terms. This usually involves the addition of prior 
information on the other cosmological parameters. For the 
work presented here, this means in general to add a prior 
from the forthcoming Planck CMB surveyor and this probe 
will be discussed in the next section. This means we just 
have to multiply the likelihood p with the prior to obtain 
the posterior, or in terms of the Fisher information matrix 



+ F° 



(4) 



The remaining procedure is to integrate the posterior over 
the nuisance parameters and the cosmological parameters 
not involving the equation of state. In terms of the Fisher 
matrix approximation this can be achieved by inverting F, 
projecting this inverse to the space involving the equation 
of state parameters Wi and inverting again, i.e. 

-I -1 



(F-) 



(5) 



where the index {i,j} runs over the indices of the 9k corre- 
sponding to the w bin parameters. 

The matrix F is the estimator of the correlation matrix 
we want to decorrelate. We will hence calculate the eigen- 
values Xi and the eigenvectors ei{z), where the vector has 
the dimension of the number of bins N. In the basis of the 
eigenvectors 6^(2) the Fisher matrix is diagonal and hence 
decorrelated. The eigenvectors represent the principal com- 
ponents of the corresponding survey. We can now write the 
underlying w{z) in terms of the eigenvectors 

N 

w{z) = '^aiei{z), (6) 



of the eigenvalues Aa,; 
given by 



a; 



-1/2 



Hence the error on w is 



dw{zf = ^ Aa'^Biizf = ^ 
1=1 1=1 



!i(£)! 



(7) 



where ei{zy = {ei^^ea, ■■■, e^jv)^ is a vector 

Since we have now a hierarchy of modes, with increas- 
ing error bars, we can in fact use a criteria similar to Oc- 
cam's razor to decide how ma ny modes are constrained sig - 
nificantly by the experiment. iHuterer fc StarkmanI (|2003 ) 
used a risk factor to decide this. This includes the increas- 
ing variance with the number of modes, but also a decreas- 
ing bias factor. In a modern Bayesian approach this prob- 
lem can be addre ssed with the Bayesian informa tion crite- 
rion or evidence (jSaini et al.l 120041 : iLiddld IJOmI ') . We wiU 
come back to this in the discussion section of the paper. In 
the meantime we will concentrate on studying the structure 
of the eigenmodes and thei r associated eigenvalues. While 
lAlbrecht fc BernsteinI (|2007r ) concentrated on the area of the 
posterior probability for a given number of modes to assign 
a figure of merit to an experiment, we will in addition study 
the redshift dependence of the eigenvectors. This will allow 
to shed light on the redshift sensitivity of a given survey. In 
order to study this dependence and accuracy we will plot in 
the forthcoming sections the quantity 



(8) 



(t>i{z) = A VAiei(z) 



where the amplitude of this expression is representative to 
the accuracy of the mode and ei{z) is encoding the redshift 
sensitivity. Note that the prefactor A appears in Eqn. ^ in 
order to make this quantity quasi-independent to the num- 
ber of bins. 

We will now proceed to forecast the parameter con- 
straints from the Planck surveyor, which we subsequently 
will use as prior information for the other surveys. 



3 COSMIC MICROWAVE BACKGROUND 

In this section we describe our principal component anal- 
ysis for the case of the CMB as a cosmological probe. 
The Fisher matrix for CMB power spectrum is g iven by 
(jZaldarriaga fc Seliaklll997l : IZaldarriaga et al.lll997l ) 



F"- 



EE^cov-^^^ 



I X,Y' 



dOi 



(9) 



where the variables X and Y represent different cross and 
auto correlations in the CMB power spectrum, i.e. Cx,i rep- 
resent the TT, EE or TE power spectra where T repre- 
sents temperature anisotropy and E represents E-mode po- 
larization anisotropy. The covariance b etween auto and cross 
power spectra is given by COVxv (|Zaldarriaga fc SeliakI 
1 19971 : IZaldarriaga et al.lll997l ) . 

The covariance COVxy will depend on the experi- 
mental noise of a given CMB experiment. In this paper, 
we choose the Planck experiment as a baseline to define 
the noise attributes that will determine COVxy. We re- 
fer the reader to the Planck blue book for these values 
(|ESA-SCI(2005Ull2005l ). The Planck mission is designed to 
measure the CMB anisotropics in nine frequencies over the 



Future Probes of Dark Energy 5 



Parameter 


Fiducial Value 


n^ 


0.3 


fib 


0.04 


0-8 


0.78 


Ho 


72 


Us 


1.0 


T 


0.09 


u,(z) 


-0.9 



Table 1. The fiducial cosmological model used for the the fore- 
casting analysis in this paper. 



entire sky. We assume here that the different frequencies will 
provide enough information so that foregrounds can be re- 
moved over as much as 80% of the sky in order to perform 
cosmological analysis over this area. We assume only one 
science frequency is used after foreground cleaning, e.g. the 
143 GHz channel. 

Usually one can use all the information including low 
multipoles to determine cosmological parameters from the 
CMB. However, in this paper we are concentrating on the 
role of the dark energy equations of state. In order to treat 
low CMB multipoles correctly one requires to calculate the 
integrated Sachs- Wolfe (ISW) effect. For a consistent treat- 
ment of the ISW it is necessary t o include dark energy per- 
turbations IjWeller fc LewisI 120031 ). There is a singularity in 
the perturbation equations (not necessarily essent ial) , when 
the equation of state c rosses the line w — —1 (jHul l2005l : 
ICaldwell fc Doranll2005l ). which of course for arbitrary bins 
can happen. In addition there arises a problem for the bin- 
ning approach in context of the ISW. Since the perturba- 
tions in the dark energy component depend on the time 
derivative of the equation of state, any step binnin g in w 
involves singular derivatives l| Weller fc LewisI |2003| ) . This, 
of course, could be resolved by imposing continuous bins 
ICrittenden fc Pogosianll2005l ). In light of these difficulties 
we choose not to include CMB multipoles below £min = 20. 
The drawback of this approach is that the constraints from 
the CMB on the optical depth are poor, but we prefer to 
have a conservative estimate rather than treat this pertur- 
bations incorrectly. 

If we ignore foreground noise and the late time ISW ef- 
fect, Eqn. (|9]) forecasts the errors on parameter set 9i given 
by the Planck experiment. For the purposes of this paper we 
will treat Eqn. ([SJ as a prior for other dark energy probes. 
Note that in addition to the cosmological parameters intro- 
duced in section (2] the CMB analysis requires the optical 
depth T due to reionization of the Universe as a free param- 
eter. 

The fiducial values of the cosmological parameters we 
use for the analysis are listed in Table ([1} . In the remainder 
of the paper we refer to the cosmological parameters ex- 
cluding the w-bin parameters as standard cosmological pa- 
rameters. In this p aper we use a modified version of CAMB 
l|Lewis et al.l 120001 ) to calculate CMB power spectrum; we 
modified the code by allowing for binning in w as defined in 
Eqn. IT}. In the case of the CMB alone we choose Zmax = 3. 
This has been chosen at a relatively low redshift compared 
to the CMB redshift because for a general equation of state 
w^z) < it is mainly the low redshift behaviour that has 
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Figure 1. We plot <j>i{z) = N\/\{\ei\ for the Planck experiment. 
The solid line indicates the first best estimated principal compo- 
nent (PC). The dotted line shows the second best estimated PC 
and the dashed line indicates the third one. The faint dotted lines 
show the remaining PCs. The upper panel represents the PCs 
with fixed cosmological parameters, while the lower panel shows 
the PCs after cosmological parameters have been marginalized. 



the h ighest impact on CMB observables IjDeep Saini et al.l 
l2003l ). The choice of a higher Zmax > 3 will have very lit- 
tle impact on the principal components that we study. For 
the convenience of joint analysis on different probes, we set 
the redshift width of each bin as a constant Az = 0.05. We 
emphasize that our results do not change if we decrease the 
size of the bins. With this set up it is simple to perform a 
joint analysis between two experiments. Even though two 
experiments may have two different maximum redshifts and 
hence a different total number of bins N, in the overlapping 
redshift range the wi axe the same and it is straight forward 
to add their Fisher matrices and combine them to produce 
principal components for joint experiments. In order to do 
this, we simply add corresponding columns in the Fisher 
Matrices. For parameters, which are only relevant for one 
experiment but not for a second experiment, we insert ze- 
ros in the Fisher matrix given dOj/dOi = for the relevant 
observable Oj and 9i for the second experiment. 

In Fig.[l]we plot (ji-ii^) a-s given in Eqn. (|8]) for the Planck 
experiment. The upper panel shows 0i(z) when cosmolog- 
ical parameters are fixed. In this panel, the first and sec- 
ond components are three orders of magnitude higher than 
the rest of (j)i. No sign transition occurs on the first mode, 
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while on the second mode <?i2 changes sign around z = 0.25. 
This is indicated by the kink in the logarithmic plot. In the 
lower panel we plot <j>i{z) after the standard cosmological 
parameters have been marginalized. Compared with the up- 
per panel, all (f>i are one order of magnitude smaller. This 
is due to the degeneracy between ■w{z) and the cosmological 
parameters, in particular fim. 

To visualize the impact of the principal component on 
the observable, we plot in Fig. [5] the change of Cx.i when 
wi^z) changes along the direction of the principal compo- 
nents, i.e, we construct w{z) as 



w{z) = Wfid(z) + AQiei(2:) 



(10) 



where Aoi is the error on a\. For a data set with 60 param- 
eters, Aai = 8.2A/Ai~ give the 1-a boundary along the i-th 
eigenmode direction from the center. Note that the prefac- 
tor in this relation is there to represent a l-cr errorbar in a 
60-dimensional parameter space. 

Figures [5] (a), (b) and (c) presents the changes of Ctt, 
Cbe and Cte, respectively. The solid and dotted line repre- 
sent changes in the first and second eigenmodes, respectively. 
In the top plot of each panel, we show how w deviates from 
the fiducial value (faint dotted). For i = 1 (solid), w only 
deviates slightly from wnd at z < 0.5; while for i = 2 (dot- 
ted), w descends from —0.6 at 2: = 0. to —1.1 around z = 0.6 
and then goes slowly back to —0.9 at higher redshift. The 
second plot from the top in each panel shows how Cx,i be- 
haves. The light area shows the observational error. For ev- 
ery spectrum, the deviation from the fiducial Ci is very small 
compared with the observational error, in fact too small to 
be seen in the graph. For clarity, we presents in the bottom 
plot in each panel the relative difference ACx,i/Cx,i to show 
how Cx,i changes. The peaks of ACx,i/Cx,i have the same 
height, which indicates that Cx,i shifts along / when we add 
the principal components. This is because ■w{z) changes the 
expansion history of the Universe and therefore the angu- 
lar positions of acoustic peaks. Ctt, Cee and Cte all have 
the same qualitative behaviour with changing eigenmodes. 
Notice in the bottom panel of (c) that the peaks become 
infinity because Cte crosses zero at the same position. 

The 'constraint' on w{z) from CMB power spectrum 
mainly come from the estimation on the angular posi- 
tions of acoustic peaks. Since the angular positions also 
depend o n ^m, there is a large degeneracy between Sim 
and tofzlte ean fc Melchiorrl 120021 '). Therefore, X^^ in the 
marginalized case is two orders of magnitude smaller than 
the one with fixed cosmological parameters. However, the 
first mode is still three orders of magnitude higher than 
the rest of (pi. The second mode drops down two orders of 
magnitude, but it still changes sign around z = 0.25. The 
lower panel is consisten t with the lowest panel of Fig.(l) in 
ICrittenden fc PogosiaJ (120051 ). 

From the discussion above, one can tell that the infor- 
mation for w from CMB power spectrum is very limited. 
However, since the Planck surveyor will be able to pin down 
the other cosmological parameters to a percentage level, we 
will choose it as a prior on cosmology when we evaluate w 
from other dark energy probes. 



4 TYPE lA SUPERNOVAE 

Type la SNe are so called standardis able candles and can be 
used to probe co smological models (|Perlmutter et al.lll999l : 
iRiess et al.l Il998l ) . Their luminosity distance-redshift rela- 
tion provides a straightforward way to measure the expan- 
sion rate of the Universe. The magnitude-redshift relation 
for Type la SNe is given by 



m = M + 5 logio dL (z) + 25, 



(11) 



where m is the apparent magnitude and M is the intrinsic 
magnitude of the Supernovae. In order to forecast a given 
SNe survey, we have to ch oose a fiducial value fo r M, and 
we assume Mfld = -19.3 (|Perlmutter et al.lll999l ). The lu- 
minosity distance dj^ in a fiat Universe is given by 



di^iz) = c(l + 2:) 



r dz' 

Jo H{z' 



(12) 



where c is the speed of light and H{z) is the Hubble param- 
eter. For a fiat universe, H{z) is defined as 



H{z) = Ho 



nn,{i + zf + {i 



ilm) exp 3 / ——^ — -—dz 
\ Jo 1 + z' 



(13) 

The Fisher ma trix for the Supernovae survey is given 
by l|Tegmark et al.i ri998l 



F ■ — 



Jo '^m 



dnia dnio 



■n{z)dz, 



(14) 



where n{z) is the density distribution of the Supernovae sat- 
isfying 



n(z)Az — N, 



(15) 



where A'' is the total number of SNe la in the survey, Zmax de- 
notes the survey depth and am is the error on the magnitude 
m. We assume that am = 0.15 for all Type la SNe surveys, 
which is a simplistic, but common assumption neglecting, 
for example, measurement errors and cadence differences. 

We start our discussion with the SN AP survey 
||SNAP Collaboration! [2OO5I : [Albert et al.ll2005^ . SNAP is a 
proposed space mission designed to measure the light curves 
and spectra of Type lA SNe and the spectra of their host 
galaxies to estimate their redshifts. It is estimated that up 
to 10, 000 SNe la could be found out between z = 0.3 and 
1.7, 2, 000 out of which will have well-measured light curves 
(|SNAP Collaboration! I2OO5I : [Albert et al.ll2005l ). Note that 
the value of Zmax in Eqn. ^ will be different depending on 
different survey parameters and for the SNAP SNe survey 

IS S^max — i . ( . 

For a fiat Universe the only standard cosmological pa- 
rameters relevant for Supernovae as cosmological probes are 
the total matter density f2m and Hq. The intrinsic magni- 
tude M can be combined with Ho and be treated as a 'nui- 
sance' parameter and marginalized over. However, we treat 
both of them separated given that Ho can be constrained to 
a certain extent by the Planck surveyor. 

We conservatively assume that we have 50 calibration 
SNe la between z = — 0.3. In addition, we assume that the 
rest of SNe la are uniformly distributed in redshift. 

The evolution of the eigenmodes in z, therefore, will 
mainly depend on the derivatives of the apparent magni- 
tude m{z) with respect to Wi. Fig. Q shows all (jf>i for this 
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Figure 2. The changes of Cx i as we perturb w around the fiducial model along the eigenmode directions by l-tr with standard 
cosmological parameters. In (a), (b) and (c) we show the changes of Ctt, C'ee and CtEi respectively. The solid and dotted line show the 
first and second eigenmode, respectively. The faint dotted line is the fiducial model. In the top plot of each panel, we show w{z) given 
the change along the different principal components. In the middle plots the light area shows the observational error on each spectrum 
per multipole /. In the bottom we show <5Cx.i/Cx,i- 
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Figure 3. The (pi for the SNAP SNe la survey. The solid line 
shows </>i(z). The dotted line and dash line show <j>2{z) and <f>3{z), 
respectively. The faint (green) dotted lines show the remaining 
modes. The upper panel is for fixed cosmological parameters, 
while the lower panel is after cosmological parameters have been 
marginalized including the prior from Planck. 



survey. There are 34 eigenmodes in Fig. Q. As in the CMB 
case, we highlight the first three eigenmodes. The (black) 
solid line, the (red) dotted line and the (blue) dashed line 
represent 4>i, 4'2 and 03, respectively. The faint (green) dot- 
ted lines show all higher order eigenmodes. The upper panel 
shows (f)i maintaining fim and Ho fixed, while the lower panel 
represents the eigenvectors after marginalizing over Jim and 
Ho including the Planck prior. The errors on the compo- 
nents drop down linearly from the best constrained ones to 
the worst ones. By marginalizing over fim and Hq, the am- 
plitudes of 4'i become slightly smaller, w hich is consistent 
with the results of Hutcrcr fc StarkmanI l|2003r ). We notice 
that the redshift dependence of the eigenmodes has changed 
with marginalization; the eigenmodes acquire zeros as we 
marginalize over ilm and Hq, indicating the impact of pa- 
rameter degeneracies. 

Fig-Hlshows how the apparent magnitude m{z) changes 
when w changes along the direction of the eigenmodes. We 
take the first three components as examples. In the upper 
panel, we plot w — uifld + Aciiei(z) where Aai = 6.1A|~ 
for the 34 free parameters. The middle panel shows the mag- 
nitude m{z). For clarity, we also show in the lower panel the 
absolute change \5m\ relative to (Jm = 0.15 for each case. 
The solid line is when we vary the first eigenmode, while 
the dotted and dashed lines represent variations of the sec- 
ond and third eigenmodes, respectively. As we expect, the 
error on m due to these variations is less than (Jm = 0.15. 

We can now proceed to compare the (f>i from SNAP, 
PS4 and DES SNe la surveys. PS4 will have a medium 
survey in the grizY bands over 1,200 deg^ and an ul- 
tra deep survey in the same bands over 28 deg^, which 
is two magnitudes deeper than the medium survey. These 
will allow potentially the detection of 5, 000 Type la SNe 
and 1,000 in early type hostqj. We note that DES and 
PAN-Starrs will not have the capability of spectroscopic 
follow-up. Therefore follow-up will have to be done with 
other instruments or photometric redshifts will have to be 



see at: pan-starrs. ifa.hawali.edu/public/science/active.html 
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Figure 4. The change of ra(z) when w changes along the di- 
rection of the principal components. The upper panel shows 
w = Wfld + 6.1A. ' ei(z) with i = 1,2,3. The second panel 
from the top shows the m{z) for each case. In the lower panels we 
show the absolute change \5m\ relative to a^i = 0.15 correspond- 
ing to each case. The solid, dotted and dashed lines represents 
i = 1, 2, 3, respectively. 
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Figure 5. 4>i (upper panel) and 02 (lower panel) for different 
SNe la surveys. The solid line, dotted line and dash line indicate 
SNAP, DES and Pan-Starr4, respectively. 



Table 2. Experiments' Parameters For Type la SNe surveys. 
Note that we put 50 SNe la between z = and 0.3. The rest is 
uniformly distributed in redshift. 



used. We do not consider these caveats in this analysis 
and assume real spectra are available for redshift determi- 
nation. DES has a dedicated Supernovae program in con- 
junction wi th a spectroscopic follow-up for a sub-sample t o 
test puritv (|The Dark Energy Survey CollaborationI 120051 ). 
Table (|4]) lists the corresponding survey parameters and the 
number of bins that we use for each survey. 

Fig. ([5]) shows 01 and (j>2 for each of these surveys. The 
solid , dotted and dash hue indicate SNAP, DES and PS4, 
respectively. Given the state of current Supernovae surveys, 
DES, SNAP and PS4 will improve on the current situation 
by providing a very long exposure time on red filters, hence 
being able to find Supernovae at higher redshifts. We have 
not considered the PSl Supernovae survey as it will not be 
considerably deeper in the red hence not providing much 
larger redshift coverage for the Supernovae detected, com- 
pared to current surveys. Due to the different distributions 
of the SNe and different minimum and maximum redshift, 
the eigenmodes peak in different locations. The relative am- 
plitudes clearly show that surveys with more S Ne la give 
better constraints, as expected and also found bv lKim et al.l 
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Figure 6. Comparison of (pi between different SNe la redshift 
distribution n{z) for SNAP SNe la survey. The solid line repre- 
sents ipi by assuming constant n{z) an d the dotted line shows (pi 
with the distribution from Table(l) in lKim et al] |200l). 



(|2004 ). On the other hand, DES and SNAP have a similar 
number of SNe la, but the constraint from DES is weaker 
because of its lower maximum redshift and its higher mini- 
mum redshift. 

In reality, the SNe la distribution n{z) is unlikely to be 
constant; instead, n{z) is a result of the survey strategy. In 



order to find out how the </'i(z) depend on n{z), we adopt a 
fiduci al SNAP distribution given by Table(l) in iKim et alj 
l|2004f l. Fig. ((6| shows the (j>i with different assumptions on 
n{z). The solid line in this plot represents (pi by assuming 
constant n( z), while the dotte d line shows (pi with the distri- 
bution from lKim et alJ (|200J). Beyond z = 1, the difference 
between these two cases is small. The amplitude of (pi for the 
more realistic distribution is lower than the constant n{z) 
case due to a lower number of Type la SNe at low redshift. 
The peak of (pi for the more realistic distribution is located 
at higher redshift. 
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is given by 



-Tij — 



V^ v^ 9Pa,i r„ , 

EE^^pv, 



AB] -1 OPb 



(19) 



; A.B 



with a ^ (3. Setting A = (q, /3) and B = ( a , f3 ), the ele - 
ments of the covariance Covi are given by (JMa et al.ll2006n 



Covf 



(■PA,l,ft,l) 



^/^ ^'^^la r-('i ~r ^(^ (^'\^ia ^' \ 



{21 + l)/.ky L (-•"') ("■'3') ^ ^(-./3')-^(/3.«')J 

(20) 



in which P is 



5 WEAK LENSING TOMOGRAPHY 

Weak lensing tomography is a recent technique which as- 
sumes that light from distant galaxies travels to us in 
geodesies which are perturbed by the network of dark mat- 
ter present in the Universe. These perturbations have the 
net result of changing the ellipticity of galaxies by roughly a 
per cent. By measuring the accurate shapes of galaxies and 
removing instrumental effects this cosmological effect is a 
nice probe of the growth of structure in our Universe IjHul 
I1999I ). The techniques to extract the shear signal by looking 
at the ellipticities are vast and different for space based and 
groun d based surveys IjHevmans et. akluOOq : iMassev et. al.l 
I2OO7I ). Different surveys will differ mainly by the quality of 
their image, which will affect the quality of the shear mea- 
sured, the quality of the photometric redshifts, which will 
be determined by the integration in each of the photometric 
bands and which will determine the ability to extract red- 
shift information via tomography, and the total number of 
galaxies used for weak lensing which will be determined by a 
combination of the depth and the mean seeing of the survey. 
Cosmology is extracted by the means of the weak 
lensing power spectrum. The weak lensing power spec- 
trum is the weighted two-dimensional projection of 
the three-dimensional matter power spectrum, giv en by 
iBartelmann fc Schneide3l200ll : kefregier et"aLll2004l ) 



-P{a,/3)(0 — -P{a,/3),1 + <5(a,,/3) 



(21) 



_9_ 
16 



(^)*"'"/ 



W^WpP{l/d,,z)^dz, (16) 



where A represents the pair of tomographic slices (a,/3). dc 
is the comoving angular diameter distance given by 



dc 



dL 



(1 + ^)' 



(17) 



with dL from Eqn. (|12|l . The upper limit of the integral 
Zh represents the redshift of the survey depth. The window 
function Wa{z) is given by 



Wc.{z) = - 

a 



dc(z') Oz 



Az , 



(18) 



where the weighting function W represents the strength of 
the lensing signal for a source at a given redshift and a dis- 
tribution of sources na{z) in each tomographic slice a. Note 
that to distinguish the slice index form the w-hin index, we 
use Greek letters to present the slice index in this section. 
P{l/dc,z) is the matter power spectrum at redshift z with 
l/dc representing the wavenumber k for a given harmonic I. 
The Fisher Matrix for the weak lensing power spectrum 



which is a combination of the power spectrum and the noise, 
/sky is the fraction of the sky coverage of the survey, a-, is 
the shear variance due to the intrinsic ellipticity including 
other measurement errors, rig is the average surface density 
of galaxies. Ean. (|20|) holds under the Gaussian assumption 
on the convergence field. With I SC 3000, this approxima- 
tion works verv well llwhite fc Hull200ol : ICoorav fc Hull200ll : 
iHuterer fc Taka da 2005). Therefore, in this paper, we take 
/ ^ 3000. The lowest I is limited by the sky coverage of the 
survey. Also since the cosmic variance dominates at very 
low I, low I harmonics do not contribute much to the Fisher 
Matrix, we take the lower limit of I — 10. 

One of the difficulties that we are facing in weak lens- 
ing tomography is to determine the matter power spectrum 
P{k, z) with the w-binning parameterization. Dark energy 
modifies the growth of the matter perturbation and the 
shape of the matter power spectrum. The growth factor D{z) 
is analytically given by 

d" + ^-[1- w(a){l - Q^ia))]D' - ~n^D = 0, (22) 
2 a 2 a-' 

with the prime representing the derivative with respect to 
the scale factor a. However, since future weak lensing sur- 
veys extend to the nonlinear scale, the dependence of the 
shape of the nonlinear matter power spectrum on w has to 
be considered whe n calculati ng the lensing power spectrum. 
iPeacock fc Doddsl ((1996) and lSmith et al.[ (|2003|) developed 
a set of fitting formulae for calculating the nonlinear matter 
power spectrum with a cosmological constant assumption. 
Though work has also been done to correct the fitting formu- 
lae for constant w 7^ -1 (|Ma et al.lll999l : iMcDonald et all 
2006), this area is still limited to a non-evolving equation of 
state on a small cosmological parameter space. It requires 
a lot of effort calibrating the matter power spectrum con- 
sidering evolving w by running low resolution N-body sim- 
ulations. However, since the aim of this paper is to demon- 
strate the behaviour of the eigenmodes of w with differ- 
ent future weak lensing surveys, we simply generalize the 
method given bv lSmith et al.l (|2003r ) for ly-bin parameteri- 
zation. We emphasize that in real data, either N-body sim- 
ulation with evolving w must be run or the non-linear P{k) 
must b e predicted with a renormalis able perturbation theory 
model ([Crocce fc Scoccimarrdl2006l ) or run by a training set 
i Habib et al.ll2007l ). Here we still use the fitting formulae in 
ISmith et al.l (|2003J) but use the transfer function ouput from 
modified CAMB code (see SeclS]). 
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The SNAP weak lensing survey is one of the most repre- 
sentative future surveys. Its basehne is to image over 1,000 
square degree and approximately 100 resolved galaxies will 
be found in one square arcmin (SNAP Collaboration 20051 ) 
up to redshift 3. Six optical and three infrared filters are 
used to estimate the photometric redshift of the galaxies. 
In this paper, we a dopt the simulated galaxy distributions 
in iRefregier et al.l (|2004l ). We focus on the case with the 
three tomography slices since the constraint improves in- 
significantly with more slices. We set Zmax = 3 since SNAP 
can find galaxies up to redshift 3. Fig. ([7]) shows 60 4>i{z). 
The (black) solid line represents (j)i. The (red) dotted line 
and (blue) dash line indicate 02 and <j!>3, respectively. The 
faint (green) dotted lines show all higher order modes. In 
the upper panel we fix the cosmological parameters. The 
first mode dominates at low redshifts and then decays along 
the redshift. Since the first mode drops down four orders of 
magnitude, at redshift z > 1, the second and third modes 
become dominant. The errors on the components drop down 
exponentially from the best constrained ones to those worst 
ones. In the lower panel, we show (j>i[z) after cosmological 
parameters being marginalized including Planck prior. The 
amplitude of the first mode is slightly smaller in this plot. 
However the redshift dependences of the modes change dra- 
matically. 4>i has a peak around 2: = 1. At redshift z > 2.5, 
there is less information from SNAP, so <^i is mainly driven 
by Planck Prior. (j)2 and cj)^ are one order of magnitude lower 
than (j>\ but still dominate at high redshift z > 0.8. This is 
consistent with Fig.(2) in ISimpson fc Bridle! (|2006l ) where 
they found that for cosmic shear, higher modes contributes 
significantly to weight function of Wi. 

To test how Pa,i changes with the eigenmodes, we plot 
in Fig. ([8]) Pa.i and its deviation from the fiducial value when 
we modify Wi along the direction of the eigenmodes. The 
top plot on each panel shows w = wad + Aai ei, where 

Aqi = 8i/Aj^^ for 60 free parameters. We choose the first 
three best estimated modes after marginalization as the ex- 
amples which are represented by (black) solid, (red) dotted 
and (blue) dash lines, respectively. In the second plot of the 
left panel, we show the auto correlation of the first slice, 
i.e. Pa,i with A = (1, 1). The second plot of the right panel 
shows the cross correlation between the first slice and sec- 
ond slice, i.e. Pa,i with A = (1,2). The light shaded area 
is the error around the fiducial value which is presented as 
the black solid line. For clarity, we show in the bottom the 
absolute change of |Pa,i| relative to the error bar for each 
mode. The deviation is very small compared with the error 
bar, which shows that the first three modes are well con- 
strained. However, the error in this plot only comes from a 
single shell. One also obtains constraints on w from other 
auto and cross correlations. 

In the following of this Section, we compare SNAP with 
four other weak lensing surveys: EUCLID, DES, PSl and 
PS4. SNAP and EUCLID are space based, while DES and 
PS (1 and 4) are ground based surveys. As stated before 
ground based surveys and space based surveys have differ- 
ent systematic effects which enter the weak lensing analysis. 
Image quality in space is far superior than on the ground. 
This will allow for better shape measurements in space. This 
effect is hard to include in a Fisher analysis, we would be 
able to include this effect by producing image simulation 
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Figure 7. 60 <t>i{z) for SNAP weak lensing survey. The (black) 
solid line indicates the first eigenmode. The (red) dotted line and 
(blue) dash line indicate the second and the third ones, respec- 
tively. The faint (green) dotted lines show the remaining ones. The 
upper panel is when cosmological parameters are fixed, while the 
lower panel is after cosmological parameters have been marginal- 
ized with the use of the Planck prior. 



and estimating what a^ and Tig we would have for a given 
seeing and magnitude depth. Producing image simulation is 
beyond this work, so we simply keep o-, roughly constant 
and change rig to take account for the effects of image qual- 
ity. The other main difference is the seeing. Even though 
PS4 will probe the sky to greater depths than EUCLID 
or SNAP, it will have less galaxies usable for weak lens- 
ing as the seeing lowers the surface brightness of sources 
and many faint galaxies fall below the signal to noise re- 
quired for a detection. These differences are encoded within 
Table. (O which lists the experiment parameters. To calcu- 
late the lensing power spectrum, one has to know the galaxy 
distribution. Since all these four surveys will use photomet- 
ric redshifts, one has to take into account the photometric 
redshift errors 5^. A more common method to include 5^ is 
analytical convolution of the galaxy distribution n{z) with 
the probability d istribution function of the photo-z p{z^ii\z) 
(jMa et al.l 120061 ). For simplicity, p(zph|z) is assumed to be 
Gaussian with scatter cj^{z) and bias 2;bias(-z). In this paper, 
we use mock photometric catalogues to calculate n{z) for 
each tomographic slice. These photometric catalogues have 
been produced in the following way: galaxies from GOODS 
have been cloned and their photometry has been calculated 
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Figure 8. The changes of P\ i as we perturb w around the fiducial model along the eigenvector directions. We show the change of Pa,i 
with A = (1, 1) and A = (1, 2) on the left and right panel, respectively. The (black) solid, (red) dotted and (blue) dash line represent the 
first three eigenmodes, respectively. In the top plot of each panel, we show w given by -u; = uiflj + Aoi e;. The faint (green) area shows 
the observational error on each spectrum per multipole I. we present in the bottom of each panel \5P\i/APj^i\. 



analytically. Then we have used a photometric redshift code 
to estimate the photo-z for these galaxies having set aside a 
given number of galaxies for training purposes. A more de- 
tailed discu ssion regarding the c reati on of these mocks ca n 
be found in lAbdalla.et alll|2007l ') and lBanerii et al.l (|2007h . 
Based on these catalogues, we choose five slices for EUCLID 
and PS4, three slices for DES and PSl. We also list zimax for 
each in Table. (O. 

We plot <j!>i (z) (upper panel) and 4>2 (z) (lower panel) in 
Fig. ([9}. We marginalize the Fisher matrices over the cos- 
mological parameters and lUh including a Planck prior. As 
we expect, (j>i of DES (green dashed) and PSl (blue dash- 
dotted) are about one order of magnitude lower than ifii of 
SNAP (black solid), EUCLID (red dotted) and PS4 (ma- 
genta dash-tripled-dotted). DES has a relatively small sky 
coverage and number density of galaxies. PSl covers half the 
survey, but has a relatively small number density of galax- 
ies. The redshift dependences of (piiz) are different as well. 
(j)i of DES decays along the redshift; this may be due to the 
fact that the values of the intrinsic Fisher matrix of DES is 
small so that the Planck priors have a bigger effect in the 
analysis, (pi of the rest peak around z = 0.5 but stays posi- 
tive in the whole redshift region. We notice that the redshift 
dependence of 02 do not change much for different surveys 
except for the amplitudes and the redshift where (f)2 = 0. 
We also found that with a Planck prior, (j>i depend weakly 
on how we define the tomographic slices. 

If we compare the eigenmodes between DES and PSl, 
one can also notice that PSl has better constrained eigen- 
modes than DES. Among SNAP, EUCLID and PS4, EU- 
CLID and P S4 have better cons t raints on the eigenmodes 
than SNAP. lAmara fc Refregieij (|2007l ) conclude that the 
sky coverage is the dominant parameter on improving the 
figure of merit rather than rig and Zm ■ 



Table 3. Survey parameters for weak lensing observation 



Survey 


0"7 


/sky 


"g 


2^11-1 


■2m ax 


SNAP 


0.3 


0.024 


100 


1.23 


3 


EUCLID 


0.3 


0.5 


30 


0.83 


3 


DES 


0.3 


0.12 


10 


0.67 


2 


PSl 


0.3 


0.5 


5 


0.6 


2 


PS4 


0.3 


0.5 


20 


0.83 


3 



6 GALAXY CLUSTER COUNTS 

Galaxy cluster counts probe cosmological models via the 
redshift dependence of the survey volume and the sensitivity 
of the halo mass function to the linear growth of structures. 
Forthcoming galaxy cluster counts will be performed by se- 
lecting clusters via their Sunyacv-Zcl'dovich decrement, for 
example with P lanck (|eS A-SCI(2005)1 200^) or the South 
Pole Telescope IjRuhl et al.ll2004l ) or in the x-rays with, for 
example, the extended ROentgen Survey with an Imaging 
Telescope Array (eROSITAJj. Finally future imaging and 
photometric redshift surveys will be able to identify clusters 
in th e optical wavebands ( Koes ter et al.l 120071 : IRozo et al.l 
l2007f ). For example DES and EUCLID would be able to 
identify thousands of clusters with this method. 

The redshift evolution of the number of the clusters A7V 
in the redshift interval between z — Az/2 and z + Az/2 found 
by a Galaxy Cluster survey is given by 



AiV = A^A^ 



dV 
'dzdQ 



A/ii„,{z) 



^(M,.)dM, 



(23) 



see at:|http://www.mp6-garching.mpg.de/projects.html#erosita| 
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Figure 9. ipl (upper panel) and ip2 (lower panel) from different 
weak lensing surveys. The (black) solid line, (red) dotted line and 
(green) dash line indicate SNAP, EUCLID and DES, respectively. 
The (blue) dash-dotted line and (magenta) dash-dotted-dotted- 
dotted line represent PanStarrs 1 and 4, respectively. 



where Afi is the sky coverage. dV/dzdO, is the comoving vol- 
ume element, Muin{z) is the mass limit above which clusters 
can be found at redshift z for a given set of survey parame- 
ters. dn/dM is the number density of the clusters with mass 
M and redshift z; in thi s paper we adopt the fitting function 
of dn/dM presented by I Jenkins et al.l (|200ll ) 

^(^^-) -°-3^^ A^^^ -P-|°-^^ - l«^-[^(^)-ll^'^^ 

(24) 
where pm is the matter density at present, D{z) is the 
growth factor (Jm is the linear overdensity at mass scale M 
at redshift z — 0, which is given by 



2 



W^(fc7?)A^(fc)dlnfc, 



(25) 



where R — (3/(47r)pmA/) ' . We use a spherical top hat 
window function 



W{x) = 3 



(26) 



The last term is the linear power spectrum A^(fe) = 
4nk^P{k). Following SecdSJ, we use CAMB to calculate the 
transfer function and the power spectrum P{k). 

The Fisher matrix for the cluster counts is given by 



-E 



dANc dANa 



AiVc 



(27) 



In which AN^ is the number of clusters in redshift interval 
Az at redshift Za, with z„ the maximum redshift of the 
survey. 

In order to extract cosmology from a galaxy cluster 
count survey we need to understand the relation between 
the survey parameters and the limiting mass Adiin-i{z) of 
the survey. For example for Sunyaev-Zel'dovich cluster sur- 
veys, Mii m(-g) is a function of the flux limit Sum at fre- 
quency 1/ (JBattve fc Welleiil2003l ) . In this context to estimate 
Miiin{z), one requires well-calibrated mass-temperature rela- 
tion for galaxy clusters and an understanding of the scatter 
(Lima & Hu 2004,). However in order to get an estimate of 
the ability of galaxy cluster counts to constrain the equation 
of state for different types of cluster selection, we assume a 
limiting mass, which is constant in redshift and neglect scat- 
ter. It is straight forward to include scatter in any analysis 
however it is not clear at the moment how large the intrinsic 
scatter is for different survey strategies. 

The South Pole Telescope (SPT) is a bolometric ar- 
ray observing 4,000 square degree of the southern sky in 
3-4 radio wavebands, with the main signal coming from 
the 150GHz channel. It will discover thousands of cluster 
roughly above a mass threshold of 2.4 x 10^''/i~^Mq. The 
Dark Energy Survey is observing the same area of sky as 
SPT and will provide photometric redshifts and vital weak 
lensing information for the the galaxy clusters. The observed 
redshift range goes roughly out to z = 1.5. In Fig. ([10} we 
show (j}i{z) for fixed and marginalized cosmological parame- 
ters. As before we highlight the first three components. We 
find that only the first eigenmode is dominating. This is sim- 
ilar to the SNe case, where the hierarchy of modes is almost 
linear. After we marginalize over the other cosmological pa- 
rameters, the amplitudes of the eigenmodes are one order of 
magnitude lower. Notice that the eigenmodes do not change 
significantly after the marginalization. This is because we 
impose the prior from Planck, which constrains the other 
cosmological parameters very tightly. 

In order to test the significance of the eigenmodes, we 
replace the fiducial model of w{z) with w = wgd + Aai ei 

1/2 

with Aai — 5.7A; corresponding to the 1 — a errorbars 
for 30 free parameters. In Fig. pip , we plot AN{z) for the 
fiducial model and AN{z) when w — wnd + Sai ei for the first 
three modes. The solid, dotted and dashed lines represent 
modes i — 1,2,3, respectively. In the top panel, we show 
w{z). In the second panel from the top, we show AN{z) 
together with Poisson errors around the fiducial model (faint 
dotted). For clarity, we show the absolute change |(5(AA'^)] 
relative to the error bar y/AN in the lowest panel. For some 
redshifts where z < 0.5, 5{AN) is larger than expected for 
the third eigenmode indicating that there might be redshift 
ranges where we can constrain up to three modes, although 
overall the constraint on this mode is much weaker than 
the on the first eigenmo de. This agrees with earlier findings 
(jBattve fc Wellerl 120031 ) that cluster counts will be able to 
give a good constraint on w since w can affect both the 
volume element dV/dzd^ and the growth factor D{z). 

In order to compare different surveys we show in Fig. 
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Figure 11. The change in AN{z) when w changes along the 

principal directions of the eigenmodes. The upper panel shows 

— 1/2 
w = uifld + 5.7A. ' e; with i = 1,2,3. The second panel from 

the top shows the AN{z) for each case, where the faint dotted 

line is for the fiducial model including Poisson errors, the solid 

line is for the first eigenmode, the dotted line for the second and 

dashed for the third. In the lower panel we show S{AN)/\/AN 

for each mode. 



Figure 10. (p{z) for a cluster survey with a mass limit of 
Ri 2.4 X lO^*/i~^M0, which corresponds roughly to the limit of 
the SPT survey. The (black) solid line indicates the first eigen- 
mode. The (red) dotted line and (blue) dashed line indicate the 
second and the third one, respectively. The (green) faint dotted 
lines show the remaining 27 modes. The upper panel is for fixed 
cosmological parameters, while in the lower panel wc marginal- 
ize over cosmological parameters under the assumption of prior 
information from the Planck surveyor. 



(|12p (^1(2) and 4>2{z) for three different cluster surveys 
with mass limits of Mumiz) = 2.4 x 10^'^h~'^MQ (solid), 
Miim{z) = lO"/i-'M0 (dotted), Mumiz) = 5x IO^/i-^Mq 
(dashed), all over 4,000 deg^ out to a maximum redshift of 
z = 1.5. The lower mass limit could correspond to optical 
cluster selection like DES or EUCLID, albeit any realistic 
treatment should include scatter in the mass limit for these 
surveys. Nevertheless we clearly see that a lower mass limit 
leads to much better determined modes, due to the fact the 
number of observed cluster increases dramatically for lower 
mass limits. 



7 BARYON ACOUSTIC OSCILLATIONS 

Baryon Acoustic Oscillation (BAG) arise from the fluctua- 
tion due to the sound waves in the plasma which was com- 
posed of tightly coupled CMB photons and baryons before 
recombination. After combination, photons decoupled from 
baryons; thereafter the acoustic waves were imprinted as 
fluctuation in the density fields of both photons and baryons. 
Though baryons followed the structure evolution of the un- 



derlying dark matter field, the spatial scale of the acoustic 
oscillation is well reserved as a standard ruler which can be 
used to calibrate the geographical evolution of the Universe. 
As galaxies are tracers of baryons, galaxy redshift surveys 
have been performed to be a complementary probe to con- 
strain c osmological parame ters and the properties of dark 
energy (|Percival et al.ll2007l ). 

The basic idea of BAG is to calibrate the length scale of 
the acoustic waves at different redshifts and fit this distance- 
redshift relation with cosmological models. The numerical 
details of the methodology vary according to different au- 
thors. In this paper, we use the methodology developed in 
ISeo fc Eisensteinl l|2003) . where the full matter power spec- 
trum P(k, z) has been used to measure Hubble parameter 
H{z) and the angular diameter distance da(z) given by 



d^{z) 



(l + zy 



(28) 



P(k, z) is estimated from galaxies in a redshift shell cen- 
tered at z. The fisher matrix at the shell is given by 
l|Seo fc Eisensteinll2003l ) 

dlnP(k,ii)dlnP{k.fi),^ ,, .2-Kpdkdfi 

VcB(K,fI)- 



dOi 



dOi 



2(27r)3 



(29) 
where k is the norm of the wave vector k, fi is the cosine of 
the angle between k and the line of sight. T4ff is the effective 
volume of the redshift shell given by 



Veff = 



n{v)P{k,ii) 



n(r)P(fc,^) + 1 



dr, 



(30) 
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Table 4. Survey Parameters for BAO 
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Figure 12. The first (top panel) and second (lower panel) 
weighted principal component 0i {z) and 02 (2) for three different 
mass limits of Mii^{z) = 2.4 X W^'^h~'^MQ (solid), Mii^{z) = 
IQI^/i-IMq (dotted), Mii^{z) = 5 x lO"/i-iAf0 (dashed). 



in which r is the vector along the line of sight and n(r) is 
the comoving number density of the galaxies. 

The observed power spectrum P{k, /i) in Eqn. 
given by (.Sco fc Eiscnstcin 2003) 



P{K. 



dl(z)rciH{z) ^2 
dlH{z)ref 



l+P 



/C I ~r /C|| 



xD\z)P{k) + Psi 



(31) 



in which 1. and || present the traverse and line-of-sight direc- 
tions, respectively. The index "ref" indicates the reference 
cosmology we use. In this paper, we set the reference cos- 
mology the same as the fiducial model that we assume. Pshot 
is the shot noise power, b is the bias factor and (3 is ^2^/'^- 
We treat the bias of each shell independently as an nuisance 
parameter and marginalize over them in the analysis. In re- 
ality, the fiducial value of the bias depends on galaxy types 
that one uses to recover the power spectrum. Here we as- 
sume b = 2 for all galaxy surveys expect for the SKA, where 
we assume b = 1. 

One uncertainty in Eqn. (|29|l is the upper integral 
limit of k. This has been discussed by different authors 
iBlake fc Glazebrook|[2003l : ISeo fc EisensteinI |2003| ) . In this 
paper, we choose a conservative value of fcmax = 0.2hMpc~^ 
for shells at z < 2 where most surveys reach. For WFMOS 
deep, we choose fcmax ~ 0.5hMpc~^ , as it is a much higher 



Survey redshift range /gj^y 



N. 



gal 



WFMOS(wide) 

WFMOS(deep) 

SKA 



(0.5,1.3) 

(2.3,3.3) 

(0.,2) 



0.05 2 X 10^ 

0.0075 6 X 10^ 
0.5 1.4 X 10^ 



redshift survey and will probably have non-linearity which 
will prevent us to measure P{k) for larger k. 

For spectroscopic galaxy surveys, if we ignore the cor- 
relation between the shells for very small k along the line 
of sight, the Fisher matrix of the whole survey can be ob- 
tained simply by adding together the Fisher matrices at dif- 
ferent redshift shells. However, one has to be careful with 
the photometric redshift surveys. The relative large photo-z 
error will smear the power spectrum along the line of sight 
and introduce a cross correlation between different shells. In 
this paper we ignore the cross correlation between different 
shells, but given that the photo-z damps the power spectrum 
along the line of slight, i.e. 

Pphotoz = P(fcrof±, fcrofll) exp(-fcrof||'7zC/i/(z)), (32) 

where CTz is taken from simulations (jAbdalla et al.l 120071 ) . 
We choose the width of the shells to be much larger than 
the photo-z error a^ of galaxies in that shell. This strategy 
is also a solutio n for avoiding the cross correlation between 
shells (jSeo fc E iscnstcin 2003). 

Another freedom in the analysis is how to choose the 
redshift shells. The principle is to make sure that the width 
is large enough that one can measure two or three wiggles in 
the power spectrum along the line of sight in that shell. At 
the same time, we wish to make the shell thin enough so that 
we have enough information for Wi at each bin. Therefore a 
choice of Az — 0.2 as the shell width means reasonable for 
most surveys. For WFMOS deep, we use Az > 0.2 (see 
Table.©). 

Two out of the three surveys planned by WFMOS 
("wide" and "deep") will provi de a competitive constr aint 
on cosmology and dark energy (JGlazebrook et al.ll2005h . In 
this section, we take the wide survey as an example and 
study how (/)i behaves in this case. Table. Q shows the pa- 
rameters for WFMOS survey and in Table. ([5} the locations 
of the slices. We use 26 ui-bins for ^^ax ~ 1-3. We plot <j)i{z) 
in Fig. (|13p . The (black) solid line indicates the best esti- 
mated eigenmode. The (red) dotted line and (blue) dash line 
indicate the second and the third eigenmode, respectively. 
The (green) dotted lines show the other eigenmodes. In the 
upper panel we marginalize over the biases of each of the 
shells. The first mode dominates at low redshift but decays 
slightly along z and crosses zero around z = 0.8. At redshift 
z > 0.5, the second and third modes dominate. In the lower 
panel, we show (jiii^) after marginalization over all other 
parameters including a Planck prior. In this plot the first 
eigenmode peaks around z — 0.5 where the survey starts. 
In both panels, the value of 01 is still significant, which 
is different from o t her p robes. This is also consistent with 
ISimpson fc Bridle! (|2006l ) where they found that the weight 
function of BAO shows a high sensitivity on w{z) at high 
redshift due to the location of the data at high redshift. The 
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Table 5. The location of the redshift shells in WFMOS 



Survey 
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WFMOS{wide) 0.6 0.8 1.0 1.2 
WFMOS{deep) 2.5 3.0 
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Figure 13. 26 <t>i{z) for WFMOS large survey. The (black) solid 
line indicates the first eigenmode. The (red) dotted line and (blue) 
dashed line indicate the second and the third ones, respectively. 
The (green) faint dotted lines show the remaining modes. The up- 
per panel is when cosmological parameters are held fix, while the 
lower panel is after cosmological parameters have been marginal- 
ized with the of the Planck prior. 



second and third modes are also significant at high redshift. 
The obvious discontinuity are found at the redshift where 
the shells are located. 

In Fig. (|14p . we show how P{k, fi) changes along with w. 
We take the first redshift shell in WFMOS wide survey as an 
example. The left panel shows the power spectrum perpen- 
dicular to the direction of the line of sight which corresponds 
to /i = 0, while the right panel shows the power spectrum 
along the line of sight which corresponds to ^ = 1. In all the 
panels, the (black) solid line indicates the first eigenmode. 
The (red) dotted line and the (blue) dashed line represents 
the second and third modes, respectively. The top plot on 

each panel shows w — Wfid + Aai ei, where Aai = 5.4i /A;~^ 
for TV = 26. Here we use ei shown in the lower panel of 



Fig. (|13|l . One can notice that w deviates from the fidu- 
cial values dramatically at high redshift, which is different 
from what we found from SNe la, WL and Cluster Count. 
The second plot from the top presents the P{k, fi) divided 
by the power spectrum without baryons and normalized to 
the same value at fc = 0. The (green) error bars represent 
the error a(P) on the power spectrum, which is given by 
jSeo fc Eisensteinll2003l ) 



-iP) = T^ 



2(27r)2 



VcB V k'^^kA/i 



(33) 



In order to be consistent, we also divide (j{P) by the fidu- 
cial power spectrum without baryons. In the bottom plot we 
show the absolute change \S\ relative to the error bar. As we 
expect, AP is smaller than the error of the power spectrum. 
One can also notice that the amplitude of the power spec- 
trums changes significantly compared with the scale of the 
wiggle. 

In Fig. (|15p . we show (l}i{z) and <j)^{z) from different 
redshift surveys marginalized with Plank prior. There are 
seven surveys shown in this plot. The (black) solid line in- 
dicates WFMOS wide; the (red) dotted line represents WF- 
MOS deep. The (green) dash line represents SKA. Note that 
we have not analysed EUCLID BAOs in this paper, which 
is a strong component of the EUCLID program. Because of 
the large cosmic volume, the height of (f>i for the SKA is 
one magnitude larger than for other surveys. For WFMOS 
deep, the data comes from redshift 2.3 < z < 3.3; there- 
fore (pi dominates at high redshift. However, (pi is significant 
at very low redshift as well. Both (j>i of SKA and WFMOS 
deep cross zero around z = 1, which is different from (j>i{z) 
of WFMOS wide. Besides WFMOS and SKA, we also show 
the eigenmodes from photo-z redshift surveys. The (blue) 
dash dotted line shows DES, the (cyan) thin long dash line 
shows PSl and the thick one represents PS4. These (l)i{z)'s 
are very similar to (t>i{z) from Planck prior. This is because 
of the damping factor on the power spectrum of the line 
of sight due to the photo-z error. The values of the Fisher 
matrices are relatively small compared with Planck prior; 
therefore the Planck prior has a larger effect on (j}i{z). 



8 JOINT PRINCIPAL COMPONENTS 

We have performed PCA analysis on a few representative 
future surveys. In the following, we will compare the eigen- 
modes for different stages and discuss the joint principal 
components from each stage. All the probes discussed in this 
section are marginalized over the other parameters including 
the Planck prior. 

In Fig. p6|l . we show 4>i{z) (upper panel) and 4>2{z) 
(lower panel) for the surveys that we have analysed for stage 
HI . The (black) solid line indicates SNe la surveys with the 
thin and thick lines representing DES and PS4, respectively. 
The (red) dotted and dash lines represents WL from DES 
and PS4, respectively. The (green) dash-dotted line shows 
the cluster count result. The remaining lines show the re- 
sult for BAO surveys; the (blue) long dashed line is for PS4 
and the (magenta) thick long dash line represents DES. The 
(blue) dash triple-dotted lines show WFMOS with the thin 
and thick lines representing deep and wide, respectively. The 
amplitudes of (j>i{z) for this stage are roughly between ten 
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Figure 14. The changes of P{k, fj.) as we perturb w around the fiducial model along the eigenmodes direction after marginalization over 
cosniological parameters including Planck prior. The left panel shows the power spectrum perpendicular to the direction of the line of 
sight which corresponds to ^ = 0, while the right panel shows the power spectrum along the line of sight which corresponds to /^ = 1. 
The (black) solid, (red) dotted, (green) dashed lines represent the first, second and third eigenmodes, respectively. In the top plot of each 
panel, we show w. The second plot from the top shows P{k,fi) divided by the power spectrum without baryons. The light area shows 
the observational error on the spectrum under the same normalization, we presents in the bottom of each panel |5P/AP|. 



and one hundred, (^i(z) of PS4 for WL dominates at red- 
shift z < 1.2; while (j>i{z) of WFMOS deep for BAO is dom- 
inant for higher redshifts. The redshift dependence of 4>i{z) 
can be classified into three types. First, (j}i{z) behaves very 
similar to the first mode of the Planck prior, which is rep- 
resented by CC and BAO from DES and PS4. This is due 
to the relatively strong Planck prior for these cases. Sec- 
ond, (j}i{z) crosses zero at low redshift. (fii{z) of DES and 
PS4 for SNe la and WFMOS deep for BAO show this fea- 
ture. However for WFMOS deep this feature is already at 
intermediate redshifts around z « 1. The third type of be- 
haviour is encountered by PS4 for WL and WFMOS wide 
for BAO. 01 (z) stays positive and peaks at the median red- 
shift. For most of the probes, (j)i is significant at low redshift 
and then decays afterward, while WFMOS deep has also a 
very significant contribution at high redshift above z — 1.5. 
In the lower panel, we notice that the amplitude of the sec- 
ond mode 02(2) is already one order of magnitude lower 
than the first mode. The dominant contribution of (j)2{z) is 
still at low redshift, with the exception of WFMOS deep. 
Due to the dominance of the first mode 01 (z) the main red- 
shift dependence of a given survey is encoded in this mode 
l|Huterer fc Starkmanl^OO^ ). 

Fig. p7[> shows the two leading modes for the stage IV 
surveys in our analysis. The (black) solid line is for type 
la SNe from SNAP. The (red) dotted line shows WL from 
SNAP and the (green) dash line shows WL from EUCLID. 
The (blue) dash dotted line shows BAO from SKA. The 
amplitudes of (piiz) for this stage are about one order of 
magnitude higher than (t>i{z) in stage III; SKA for BAO 
dominates in its redshift range z < 2. (pi of SNAP for SNe 
la and BAO for SKA have a change in sign for the most 



dominant mode, while <pi{z) of both SNAP and EUCLID 
for WL have a mode which stays positive throughout. The 
amplitude of 02 (z) is about one order of magnitude lower, 
hence of the same significance as the primary mode for the 
stage III surveys. 



From the four probes that we discussed, one can find 
that for most surveys the best constrained eigenmodes dom- 
inate at low redshift. The exception are surveys which target 
objects at high redshifts as for example WFMOS deep. 0i 
from WFMOS deep stands out at high redshift in the upper 
panel of Fig. p6|l . In the following, we discuss how the eigen- 
modes behave if we combine different probes. We perform 
the joint analysis based on the two stages. For stage III, we 
combine PS4(SNe la), PS4(WL) and WFMOS deep (BAO). 
For stage IV, we combine SNAP(SNe la), EUCLID(WL) 
and SKA (BAO). We also include the Planck prior for both 
cases. In Fig. (|18|) we show the joint eigenmodes with the 
upper panel and the lower panel representing the Stage III 
and IV, respectively. We find that 0i (z) is significant within 
the whole redshift range expect around z — 1.2 where 0i(z) 
changes sign. It is interesting to note that higher order 
modes and namely the second mode fill this gap in red- 
shift. In addition boosting 4>i{z) at high redshift, WFMOS 
deep (BAO) is complementary with PS4(SNe la), PS4(WL) 
at high redshifts. For Stage IV we find that the equation of 
state will be constrained significantly at least out to redshift 
z = 1.5. 
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Figure 15. <^i(z) (upper panel) and cj>2{z) (lower panel) from 
different BAO experiments. The (black) solid line indicates WF- 
MOS wide; the (red) dotted line represents WFMOS deep. The 
(green) dashed line represents SKA. The (blue) dash dotted line 
shows DES, the (cyan) thin long dashed line shows PSl and the 
thick one represents PS4. 



9 W RECONSTRUCTION 

So far, we have compared different surveys by concentrat- 
ing on the behaviour of 4>\{z) and <j}^{z), which allows us to 
explore the redshift sensitivity on w{z) for each probe. To 
reconstruct w{z), however, one would like to include higher 
order eigenmodes to get less bias in the reconstructed equa- 
tion of state. However, if we use higher order eigenmodes 
the error on the reconstructed w is dominated by the errors 
from these orders. In order to decide how many eigenmodes 
to use for reconstruction we require the bias and the vari- 
ance to be low. Henc e, an alternative to th e standard figure 
of merit of a survey (jAlbrecht et al.l 120061 ) is to count how 
many eigenmodes can be well constrained under a certain 
criteria. 

In Fig. (|19|l . we plot the first 15 eigenvalues weighted 
by the number of bins, for the different surveys that we used 
in the joint analysis in the previous section. All the surveys 
are marginalized over other parameters including Planck pri- 
ors. The (black) solid line shows the Type la SNe surveys 
with the filled and unfilled circles indicating PS4 and SNAP, 
respectively. The (red) dotted lines represents the WL sur- 
veys, with the filled and unfilled stars indicating PS4 and 
EUCLID, respectively. The (green) dashed lines represents 
the BAO surveys with the filled and unfilled triangles indi- 



Figure 16. (pliz) (upper panel) and (t>2(z) (lower panel) for fu- 
ture surveys in stage III. The (black) solid line indicates SNe la 
surveys with the thin and thick lines representing DES and PS4, 
respectively. The (red) dotted and dash lines represents WL from 
DES and PS4, respectively. The (green) dash-dotted line shows 
Cluster count. The rest lines show BAO; the (blue) long dash line 
shows PS4 and the (magenta) thicker long dash line represents 
DES. The (blue) dash dotted-dotted-dotted lines show WFMOS 
with the thin and thick lines representing deep and wide, respec- 
tively. 



eating WFMOS deep and SKA, respectively. We also show 
the joint analysis with the (blue) dotted-dash lines; the filled 
and unfilled squares indicating stage III and IV, respectively. 
In order to explore the problem of the w-reconstruction 
quantitatively let us assume that one only uses M eigen- 
modes. The reconstructed Wrcc{z) is then given by 



:,[z) = w{z) + ^aiei{z), 



(34) 



where it; is the value around we wish to reconstruct the 
equation of state. The oii are the best fit coefficients in 
the eige nmode basis |ei(z)|. The reconstruction strategy 
given bv lHuterer fc StarkmanI l|2003l ) is equivalent to setting 
wiz) = and the expected a\ is then given as the pro j ection 
of ei(z) on the fiducial model. ICrittenden fc PogosianI (120051 ) 
define w = Wfld under the assumption that wad is close to 
the true physical model for any reconstruction in a real data 
situation and in this case the expected ai are zero. 

We now have to employ a statistical criterion, which 
allows us to gauge how many eigenmodes to use for the re- 
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Figure 17. (t>\{z) (upper panel) and 02 (2) (lower panel) for fu- 
ture surveys in stage IV. The (black) solid line shows SNe la 
from SNAP. The (red) dotted line shows WL from SNAP and 
the (green) dash line shows WL from EUCLID. The (blue) dash 
dotted line shows BAO from SKA. 



construction of the equation of state. There are three effects 
at work, which need to be considered. First the goodness of 
the fit, which, of course, improves with the number of used 
modes, second the degradation of the errorbars with increas- 
ing number of modes and finally the bias between the true 
underlying model and the model, w, around which we re- 
construct the equation of state. As we w ill show below the 
Bayes' fac tor (IJeffrevslllOSgl : lTrottall200a ) based on Bayesian 
evidence (|Sivialll996l ) provides exactly such a criterion. 

We will proceed by describing a Gaussian approxima- 
tion of the Bayes' factor, based on Bayesian evidence, to de- 
cide the number o f significant modes . We will follow closely 
the discussion in ISaini et al.l l|2004 ). Under Gaussian as- 
sumptions, which hold for the Fisher matrix approximation 
of the underlying likelihoods, the evidence for the data D 
given the hypotheses H is approximated by 



£ = P{D\H) « P{D\9i^,H) exp(-C) 



IF + PI 



-1/2 



(35) 

with F the covariance between the w-bin parameters, 0, in 
the basis of the eigenmodes, P the covariance of the prior 
on these parameters and 0l are the parameters of the max- 
imum likelihood. C is the term which encodes the overlap of 
the prior with the likelihood, i.e. this term is measuring the 



Figure 18. The joint eigenmodes for stage III (the upper panel) 
and stage IV (the lower panel) . The (black) solid line indicates the 
first eigenmode. The (red) dotted line and (blue) dashed line indi- 
cate the second and the third ones, respectively. The (green) faint 
dotted lines show the remaining modes. For stage III, we combine 
PS4(SNe la), PS4(WL) and WFMOS deep (BAO). For stage IV, 
we combine SNAP(SNe la), EUCLID(WL) and SKA(BAO). We 
also include Planck prior in both cases. 



bias between the prior and the likelihood and is given by 

P)-^F [Gl-Op) , (36) 



C=i(0L-6»p)^P(F 



where Op is the mean of the prior, which we choose to be 
w{z) or Oi =0 in the expansion given in Eqn. 1341 If we 
make the simplistic assumption that the prior is diagonal 
and the same in each bin we obtain 

P=tLi (37) 



Note that this prior, because of its diagonal nature, has the 
same form in the original w-bin basis and the eigenmode 
basis. Here we have to make an important point. We would 
like that the rms variance, Sw, on the mean equation state 
for the prior is independent of the number of bins. Hence we 
have to scale the error in each bin. Aw according to 



(Aw)^ = N{5wY 



(38) 



Finally we will assume that the maximum likelihood 
can be written as 

P{D\e^,H) =AAe-5-^' . (39) 

We will introduce the index M to define the evidence 
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Figure 19. The weighted eigenvalues for the surveys that wo 
used in the joint analysis. Every line represents a survey. All the 
surveys are marginalized over other parameters including Planck 
priors. The (black) solid line shows SNe la surveys with the filled 
and unfilled circles indicating PS4 and SNAP, respectively. The 
(red) dotted lines represents WL surveys with the filled and un- 
filled stars indicating PS4 and EUCLID, respectively. The (green) 
dash lines represents BAO surveys with the filled and unfilled tri- 
angles indicating WFMOS deep and SKA, respectively. Wc also 
show the joint analysis with the (blue) dotted-dash lines; the filled 
and unfilled squares indicating stage III and IV, respectively. 



for M modes. Keeping in mind that we work in the basis 
where the entries in the covariance matrix F are given by 
the eigenvalues we can write the evidence for M modes as 



£^ 



A/'exp 



exp 
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\ VAw^A, + 1 



(40) 



where the first terms measure the 'goodness-of-fit' of the 
M modes, the second term the overlap between the prior 
and the likelihood and the last term is Occam's razor by 
comparing the posterior to the prior volume. We can now 
construct the Bayes' factor jJeffrevdllQSgl : lTrottall26oi ) 
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(41) 



where Xm and Xm+i are the the log-Ukelihood values at the 
best fit points, if we fit for M or M -I- 1 modes respectively. 
Qi and /3i are the best fit expansion parameters for these 
two cases. 

We will first discuss the simple case with no bias. In this 
case we choose w(yZ) — wgd. Since we are expanding around 
the fiducial model we obtain for all expansion coefficients 
Qi = 0. Since there is no difference between ignoring modes 
or putting ai = the best fit likelihood between M and 
Af -I- 1 modes is exactly same, i.e. Xm = Xm+i- Hence both 



the best fit and the overlap term vanish for the Bayes' factor 
and only Occam's razor term remains. If we further assume 
the prior is wide compared to the likelihood we obtain for 
the Bayes' factor 
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log 



yx 



M+l 



NSwJ 



logi 



CM+l 



NSw 



(42) 

Hence, in this approximate scenario the evidence ratio is 
given by the ratio of the likelihood uncertainty on the M + 
1st mode, ctm+i, to the prior uncertainty on a single bin 
y/NSw. Although the equation of state at low redshifts is 
determined almost to the 10% level by current data, it is 
much more uncertain at higher redshifts. We hence choose 
5w = 1 for the rms uncertainty on the mean w averaged 
over the entir e redshift range. According to Jeffrey's scale 
(jjeffrevslll933) we have strong evidence if the Bayes' factor 
is 1-2 and substantial evidence if it is between 0.5 and 1. We 
choose 1 as our evidence level. If we employ this we find as 
a condition for strong evidence 



N\m > 100 



(43) 



The shaded region in Fig. [19] highlights the area where this 
condition is violated. All points above the shaded region are 
eigenmodes, which are significant. We obviously find that 
the joint analysis can determine more higher order eigen- 
modes tha n the individual surveys , which is consistent with 
Fig. (1) in ICrittenden fc PogosianI (|2005l ). This comes from 
the complementarity of the different dark energy probes. For 
most of the probes, the eigenvalues descend exponentially. 

We will now continue to analyse the full Bayes' factor 
expression in Eqn.|4l]to decide for how many eigenmodes we 
have strong evidence. We will discuss two cases: As before 
the unbiased case where we reconstruct the equation of state 
around the fiducial model and a biased case where we choose 
w{z) = —1. The latter is what will happen in a realistic 
situation, where the true underlying model is not known a 
priori, although one could imagine an iterative method to 
end up reconstructing around the true underlying model, 
but this is a posterior statement. 

To illustrate how w will take effect on the reconstruc- 
tion, we give a very simple example for the SNAP SNe la 
survey. We take the eigenmodes from the SNAP SNe la sur- 
vey with fixed cosmological parameters and use our fiducial 
cosmology with w{z) — —0.9. In Fig. (|20p we plot the re- 
constructed ■u;rcc(2;) for different w and M. The solid line is 
for A'l = 2 and the dotted and dash lines for M = 10 and 
M = 30, respectively. The thin (black) lines are for w — 
and (blue) thick lines represent w — —1.0. One notices that 
no matter how we choose 'w{z), Wiac{z) will be around the 
fiducial model at low redshift, while at high redshift, uircc(z) 
will be biased towards 'w{z). This is because SNe la surveys 
have very weak constraints on -^(2) at high redshift. 

Fig. (I20p also shows the oscillation on Wicc at high red- 
shift, which is an artifact of the PGA decomposition. This 
is due to the limited number eigenmodes with M < N that 
we choose to reconstruct ^1(2) with. For M — 2, Wroc has 
no constant behaviour with redshift. As we include more 
eigenmodes, w^cc starts to converge around —0.9 for red- 
shifts 2 < 1. The oscillation at high redshift are consistent 
with the results in lHuterer fc StarkmanI (|2003l ). in which the 
analysis indicates that it is not possible to recover w{z) at 
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Figure 20. An example on reconstructing w{z) from different w 
and M for the SNAP SNe la survey. The solid line shows M = 
2 and the dotted and dash lines show M = 10 and M = 30, 
respectively. The thin lines are for tD = and (blue) thick lines 
for w = —1.0. 



redshift z > 1 with SNAP SNe data. It is hence impor- 
tant to keep in mind the redshift evolution of the eigen- 
modes for the interpretation of the reconstructed equation 
of state. This simple example indicates that the choice of 
w will make a difference for the reconstructed w{z). If we 
compare the curves with M — 10, we find that the thin 
line starts oscillating around z — 0.2, while the (blue) thick 
line starts oscillating at around z = 0.8. The dashed lines 
which show the behaviour if we include M = 30 modes, 
both the thin line and the (blue) thick line start oscillat- 
ing at z > 1. However, the amplitude of the latter is much 
smaller indicating that our initial guess val ue w = — 1 is 
much closer to the true fiducial model. iHuterer fc StarkmanI 
(|2003l l use a mean square error (MSE) criteria to find the 
optimal number of modes, which minimize bias and vari- 
ance simultaneously. However this criteria does not include 
an O ccam's razor factor by compari son to the prior informa- 
tion. ICrittenden fc PogosianI (|2005l ) use a similar MSE, but 
include the prior variance. However their reconstruction is 
unbiased and they do not use MSE to decide on the number 
of significant modes, they use MSE over all modes to com- 
pare the efi'ectiveness of different probes. As stated above we 
will use the Bayes' factor to compare how many modes are 
significantly constraint for different pr obes. We believe that 
the m easur e provided by MSE of b oth Hutcrcr fc StarkmanI 
112003) and lCrittenden fc PogosianI ([2005) is included in our 
evidence expression and incorporates both principles. Al- 
though the Jeffrey's scale is to some extend arbitrary we 
can clearly identify for how many modes there is strong ev- 
idence and for which modes there is only marginal evidence 
for different surveys. 

In Table [6] we show the number of modes with strong 
evidence for the unbiased case and in brackets for the bi- 
ased case. Note that in general the number of modes for the 
biased case are larger, although there are exceptions. For 
each survey, we include a Planck prior on the cosmological 
parameters. For most surveys, M is around 2 — 4 except for 
SKA. With the joint analysis, M will be about 8 for stage 
HI, which indicates the complementarity of the probes; for 
stage IV, A4 is 12, which is mainly driven by SKA. 



Table 6 . The number of w-eigenmodcs with strong evidence ac- 
cording to Jeffrey's scale. The number in brackets are for the 
biased case, which is the realistic scenario. For each survey, we 
marginalize over a Planck prior on the cosmological parameters. 
Note that EUCLID has a BAO component, which we have not 
analysed in the paper. 



10 CONCLUSIONS 

In this paper we have studied future constraints on w-bins 
for four dark energy probes. Type la Supernovae, weak lens- 
ing, cluster counts and baryon acoustic oscillations. 

Instead of assuming that the equation of state parame- 
ter w can be modeled by a simple function with a few param- 
eters, we bin w given by Eqn. ((T} and treat each Wi as inde- 
pendent parameter. Throughout this article, we choose fixed 
Az = 0.05 to make sure that we can compare between differ- 
ent surveys. For each probe, we choose a few representative 
future surveys. We also use a prior from the forthcoming 
Planck CMB experiment as prior on all remaining cosmolog- 
ical parameters. We find that in future mainly weak lensing 
and BAO surveys are complementary, in a sense that weak 
lensing gives tight constrain at low redshifts but BAO al- 
lows one to push to higher redshifts. However the role of Su- 
pernovae is also complementary for low redshifts. The high 
redshift sensitivity of the BAO woul d allow one to study 
also early dark energy models (Caldwell et al.ll2003r ) in an 
efficient way. However at lower redshift weak lensing and Su- 
pernovae are the more efficient probes, with the exception 
of SKA BAO measurements. Typically forthcoming surveys 
constrain 2-3 modes. Since it is likely that there will be a 
multiple of stage HI probes it is sensible to look at the com- 
bined power of these probes and we find that in this case 
five modes would be constraint. We would advocate that the 
number of well constraint modes is actually a better figure 
of merit than the one put forward by the DETF, since it 
encompasses the ability of future surveys not just to mea- 
sure the equation of state at two redshifts, today and at 
the pivot point. We address the question of the number of 
modes in a B ayesian model selection way IjSaini et al.ll2004l : 
lLiddlell2004l ). A lthough this met hod has recently attracted 
some criticism (JEfstathioul 120081 ) it provides a well defined 
framework to establish the significant number of modes for 
a particular survey. It takes into account the goodness of 
fit, the bias and Occam's razor. However it is not only im- 
portant to establish the significant number of modes. In ad- 
dition the redshift distribution of the significant modes has 
to be taken into consideration to start talking about the 
complementarity of surveys. The complementarity can than 
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be exploited to constrain eight modes for the joint stage 
III probes and twelve for stage IV if we include SKA. One 
question we have not addressed in this paper is how useful 
surveys are to probe "dark energy" beyond the background 
evolution, for example by probing structure formation. It is 
likely that in this context future weak lensing, and to some 
extend galaxy cluster surveys, will play a pivotal role. To 
conclude we would argue that all four discussed probes have 
their merits, but weak lensing and baryon acoustic oscilla- 
tion surveys seem to bear the largest promise for revealing 
the nature of the accelerated expansion of the Universe. 
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